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Abstract 



We show by exphcit construction that for every solution of the incompressible Navier-Stokes equation 
in p + 1 dimensions, there is a uniquely associated "dual" solution of the vacuum Einstein equations 
in p + 2 dimensions. The dual geometry has an intrinsically flat timelike boundary segment whose 
extrinsic curvature is given by the stress tensor of the Navier-Stokes fluid. We consider a "near-horizon" 
limit in which Sc becomes highly accelerated. The near-horizon expansion in gravity is shown to be 
mathematically equivalent to the hydrodynamic expansion in fluid dynamics, and the Einstein equation 
reduces to the incompressible Navier-Stokes equation. For p — 2, show that the full dual geometry is 
algebraically special Petrov type II. The construction is a mathematically precise realization of suggestions 
of a holographic duality relating fluids and horizons which began with the membrane paradigm in the 
70's and resurfaced recently in studies of the AdS/CFT correspondence. 
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1 Introduction 



The Einstein equation 

G^. = 0, (1) 

and the incompressible Navier-Stokes equation 

V, - T^d'^v, + d,P + v^d,v, = 0, (2) 

have long played central roles in both mathematics and physics. The Einstein equation universally governs 
the long-distance behavior of essentially any gravitating system, while the incompressible Navier-Stokes 
equation universally governs the hydrodynamic limit of essentially any fluid. They both display a rich 
non-linear structure which has been a continual source of interesting surprises yet remains, centuries 
after their discovery, incompletely understood. It is the purpose of this paper to give a mathematically 
precise relationship between ([T]) and ([2| and their solutions, and thereby provide a hopefully useful bridge 
between the two subjects. For example cosmic censorship could be related to global existence for Navier- 
Stokes or the scale separation characterizing turbulent flows related to radial separation in a spacetime 
geometry, see e.g. [l]|2]. 

Hints of such a connection, summarized in the next section, have surfaced in various forms over the 
last three decades [3 -28 , see 26 29 -32] for reviews. In particular, excitations of a black hole horizon 



dissipate very much like those of a fluid [3}j5| [T2|[T3{|28[|33}|35] , and there has been recent discussion of 
a holographic duahty relating black holes and fluids [l][8}{l0l[l2)[l§{l7)[l9)[2l][2^ Inspired by these 
suggestions, in this paper we explicitly construct a map from solutions of the nonlinear incompressible 
Navier-Stokes equation to solutions of the nonlinear Einstein equation. A key ingredient is the imposition 
of boundary conditions which, in a sense to be defined, isolate the horizon dynamics from the rest of the 
gravitational dynamics and thereby reduce equation ([T]) to equation ([2]). 

Our basic construction is roughly as follows. We begin with the region of p + 2-dimensional Minkowski 
space inside a hypersurface Sc given by an equation of the form — — Arc- Sc is intrinsically flat 
(being the translation of an hyperbola in the t-x plane along the remaining p spatial directions) , but has 
an extrinsic curvature linked to the constant acceleration a = l/^/4r^. It asymptotes to its future horizon 
which is the null surface x = t. We then study the effect of finite perturbations of the extrinsic 
curvature of Ec while keeping the intrinsic metric flat. These generically lead, when evolved radially 
inward with the Einstein equation, to singularities on 'H'^ . The special ones which are smooth on 7^"*" 
are analyzed in the hydrodynamic "e-expansion" , which is a nonrelativistic, long- wavelength expansion 
and, importantly, keeps terms that are nonlinear in the size of the perturbation. It is found that tensor 
and scalar modes of the metric decouple in this limit and the remaining degrees of freedom are vector 
modes governed by the Navier-Stokes equation 'mp+1 dimensions. We present (equation ( 14 1 below) the 
p -\- 2-dimensional solution of the Einstein equation through third order in the hydrodynamic expansion 
parameter e. The flrst term is flat space. The second and third terms are algebraically constructed from 
the velocity field v'^ and pressure P of an incompressible fluid. The nonlinear spacetime Einstein equation 
then reduces to the nonlinear incompressible Navier-Stokes equation for the pair {v\P). 

This result is already interesting and non-trivial, but the fact that the Navier-Stokes arises when the 
geometric variables are subject to the same kind of expansion used in fluid dynamics might have been 
anticipated. A deeper connection appears when we consider an alternate expansion in which, instead of 
going to long distances, we take the acceleration of Ec to infinity. This is a near-horizon limit since it 
pushes Sc towards its horizon . We then show that, after a constant overall rescaling of the metric, the 
near-horizon expansion is mathematically identical to the hydrodynamic expansion. Hence the solutions of 
the Einstein equation (constrained by the boundary conditions of a flat metric on Ec and smoothness on 
'H'^) in this near-horizon expansion are in one-to-one correspondence with solutions of the incompressible 
Navier-Stokes equation. This then is the precise mathematical sense in which horizons are incompressible 
fluids. 
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It is possible that the uhimate origin of this relation is a deep and exact holographic duality relating 
(among other things) quantum black holes to fluids as has been suggested by string theoretic investiga- 
tions. However in this paper we have concentrated on simply establishing the mathematical relationship 
between and ^ in a manner which makes no assumptions about or reference to this tantalizing 
possibility. 

This paper is organized as follows. In section 2 we briefly describe precursors of our construction 
going back to the 70's. Section 3 briefly reviews the hydrodynamic expansion in the study of fluids, and 
the emergence of the incompressible Navier-Stokes equation in the hydrodynamic limit. In section 4 we 
specify the boundary conditions, explained roughly above, used to isolate horizon dynamics. In section 

5 we present the general solution of the nonlinear Einstein equation with these boundary conditions 
through the first three orders in the hydrodynamic expansion, and show that the first nontrivial term 
corresponds to the velocity field of an incompressible fiuid. We also discuss the geometric analog of forcing 
the fluid, argue for uniqueness, and discuss the possible formation of black hole type singularities. Section 

6 presents a simpler form of the metric and shows that, up to an overall rescaling and after an appropriate 
coordinate transformation, it depends only on the product of the leading-order acceleration of Sc and the 
hydrodynamic expansion parameter e. In section 7 we show that the geometries are, through the order 
constructed, of a special type known in four dimensions as Petrov type II. This may enable a connection of 
the present work with the large literature on algebraically special spacetimes 36 -38 . Finally in section 
8 we demonstrate, using the simplified metric of section 6, the equivalence of the hydrodynamic and 
near-horizon expansions. 



2 Relation to previous work 



The first suggestion of a relation between horizon and Navier-Stokes dynamics appears in the prescient 
thesis of Damour [3|. This work contains an expression now known as the Damour-Navier-Stokes equation 
39 governing the geometric data on any null surface. Although tantalizingly similar, it is not quite the 
Navier-Stokes equation as it has too many variables (eliminated herein by an appropriate boundary 
condition) and an extra nonlinear term. To get precisely Navier-Stokes we found it necessary to consider 
the near-null limit of a highly accelerated timelike surface. Such a surface was introduced by Price and 
Thorne [s] and coined the stretched horizon (analogous to our Sc). Although similar in spirit their limit 
is slightly different from ours. They obtain a compressible fluid with a negative bulk viscosity and an 
extra term not present in Navier-Stokes. This approach was developed into the membrane paradigm 
and is reviewed in the book [40]. Much more recently Policastro, Starinets and Son [o] made the striking 
observation, in the context of the AdS/CFT correspondence, that the dissipative behavior of a large black 
hole in AdS agrees with that of the hydrodynamics of the holographically dual CFT. This observation 



has been developed in many directions 10 12 291141 51 . Although far from obvious at first glance, these 



results from AdS/CFT are compatible with, and in some cases equivalent to, the earlier results from pure 
gravity [I0p3p8p2p3] . In the AdS/CFT context Bhattacharya, Minwalla and Wadia [Tp3] showed that 
in asymptotically AdS spacetimes at finite temperature, the asymptotic AdS boundary data is governed 
in a hydrodynamic limit by the Navier-Stokes equation. They use the Navier-Stokes data to construct 
a bulk solution of the Einstein equation with negative cosmological constant. Our dual bulk geometry 



in equation (14) is a refinement of expression (4.4) in |1 in which the cosmological constant is taken to 



zero and the boundary can be pushed to any radius - in particular to the interesting near-horizon region. 



Finally we rely heavily on our previous paper 28 which solves the linear case within the framework 
adopted herein. 



3 



3 The hydrodynamic limit and the e-expansion 



The incompressible Navier-Stokes equation has a well-known scaling symmetry which is important in the 
following and briefly reviewed here. Let the pair (w^, P) solve the incompressible Navier-Stokes equation 

d'v, = 0, drv, - rid'^v, + d,P + v^djV, = 0, (3) 

where 77 is the kinematic viscosity and i = l,...p. Now consider a family of pairs (w|,P') in which 
frequencies and wavelengths are non-relativistically dilated and amplitudes scaled down by the parameter 
e: 

vl{x\T)^tv,{ex\e\), (4) 
Pt{x\T)^e^P{ex\e\). 

It is easy to check that ^ directly implies 

d-rvl - r^d^vl + d,P' -t- v'^djvl = 0. (5) 

Hence Q generates from the original solution a family of solutions parameterized by e. 

In real fluids there are always corrections to the Navier-Stokes equation. Galilean invariance requires 
that these vanish for constant Vi. Typical corrections are for example of the form 

drv, - T]d^Vi + d,P + v^djv, + v'^v^dkdjVi + d^v, = 0. (6) 

If {vi,P) obey this equation, the rescaled quantities obey 

drvl - r^d^vl + d,P' + v^'djvl -t- {v'^'v'^ dkdjvl + = 0. (7) 

The limit e — >■ is the hydrodynamic limit. In this limit these corrections become irrelevant. Similarly 
the speed of sound goes to infinity and compressible fluids become incompressible. It is not hard to show 
that all reasonable types of corrections are scaled away, and the incompressible Navier-Stokes equation 
universally governs the hydrodynamic limit of essentially any fluid. The limit is an incredibly rich and 
interesting one because, even though the amplitudes are scaled to zero, nonlinearities survive. It is this 
hydrodynamic limit of a fluid that we will match to a near-horizon limit in gravity. 



4 Characterizing the dual geometries 

We seek a relation between the (p-|-2)-dimensional Einstein and (p-f l)-dimensional Navier-Stokes equa- 
tions. Of course, the former has a much larger solution space than the latter so only a special type of 
Einstein geometry is relevant. Roughly speaking, the relevant geometries are non-singular perturbations 
of a horizon. Let us now make this more precise. 

We consider geometries of the type depicted in Figure [T] with an outer "cutoff" boundary denoted 
Sc. The boundary hypersurface Ec is taken to be asymptotically null in both the far future and far past. 
In the Minkowskian coordinates dSp^2 — —dudv -\- dxidx^ , past null infinity I~ is the union of the null 
surfaces v — > —00 together with u —00 and Sc is the timelike hypersurface uv — —Arc with w > 0. 
Past (future) event horizons {T~L^) are defined by the boundaries of the causal future (past) of Sc- 
The dual geometries will be constructed in two a priori different expansions about Minkowski space: the 
near-horizon and the hydrodynamic e-expansion. Ultimately the two expansions will be shown to be 
equivalent. 

Initial data can be specified on the union of Ec and . We consider initial data which is asymp- 
totically Minkowskian and fiat (no incoming waves) on (or equivalently "Hr). On we generally 
demand that the intrinsic metric ^ab be fiat, 

lab^Tlab, a,b^O,...p (8) 
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Figure 1; This figure depicts tlie Einstein geometry holographically dual to a fluid. The accelerated boundary hypersurface 
Sc at radius r = is intrinsically flat but the extrinsic curvature is given by the fluid stress tensor. This extrinsic curvature 
leads to gravity waves which propagate radially inward. The leading-order condition that these waves do not cross the past 
horizon 'H~ of Sc at t = — oo or produce singularities on the future horizon 'H'^ at r = is the non-linear incompressible 
Navier-Stokes equation for the fluid. 



although wG will later consider "forcing" the system by perturbing jab- 

We wish to consider the general solution of the Einstein equations consistent with this initial data 
and smooth on In particular, so far we have not specified the extrinsic curvature Kab on Sc or 

equivalently (and more conveniently) the Brown- York stress tensor on 

Tab = 2i-fabK~Kab). (9) 

If no initial data were prescribed on Z^, any Tab on Sc consistent with the constraint equations could be 
chosen. This data could then in general be evolved radially inwards to produce a spacetime everywhere 
inside of E^- In general, such a spacetime will have gravitational flux (if not singularities) going up to 
11 = 00 (I^) as well as down to I~ . Hence we have a "shooting problem" to find those special allowed 
choices of Tab which produce a spacetime smooth on with no flux coming up from . 



We solved this problem in 28 to leading order in a double expansion in long wavelengths and weak 



fields|f] In going Rindler coordinates were used for which the leading order flat metric is 

dSp^2 = -rdr'^ + 2dTdr + dxidx\ (10) 

Ec is the accelerated surface r — Vc, is t = — oo and is r = 0. These coordinates are conve- 
nient for analyzing smoothness on It was found that the allowed choices of Tab are precisely those 
corresponding to the linearized fluid: 

r3/2T" = i)\ rl''^T^ = -■qd^'v^\ (11) 

where the (kinematic) viscosity here is given by the formula 

V = rc. (12) 



^Here we allow for incoming flux where X~ meets Sc at u = — oo, ti = 0. 
^Our normalization here agrees wi th t he conventional one for G = l/167r. 
''Our conventions here differ from 28]. 
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while Vi obeys the Unearized incompressible Navier-Stokes equation 



dy = 0, drv' - r]d^v' = 0. (13) 



If we choose any value for the viscosity other than (12), the constraint equations on are still obeyed, 
but gravitational waves are propagated down to I~ and there is a singularity at r = 0. 

In this paper we go one step further and solve the problem in certain hydrodynamic and near-horizon 
limits without making a linearized approximation, enabling us to see a direct connection between the 
nonlinear structures of the Navier-Stokes and Einstein equations. 



5 Nonlinear solution in the e-expansion 



In this section we will improve on the analysis of 28 by solving the shooting problem in the long 
wavelength e-expansion without a simultaneous linearized expansion. The general solution will be pa- 
rameterized by a solution Vi{x^,T), P{x^,t) of the full nonlinear Navier-Stokes equation with viscosity 



( 12 1 together with the parameter e. 



5.1 The solution 

Consider the metric 



rdr'^ + 2dTdr + dxidx^ 



-2(1- — ) Vidx'dT - 2 — dx'dr 



1 - 



(^2 _ ^2^ 



+ 2P)rfr2 + -^dx'dx^ 



d'^Vidx'-dT+ ... 



(14) 





2P\ 







drdr 



where Vi = Vi[x^, r) and P(x*, r) are independent of r. Here and henceforth i,j = 1, ..p indices are raised 



and lowered with Sij and we take 



(15) 



as in the hydrodynamic scaling of section 3. It follows that the first line on the right hand side of ( 14 1 is 



0{e'^) and each subsequent line is one higher order in e. The linearization of this expression in Vi agrees 



with the linearized solution studied in 28 



On the cutoff surface Ec, r — rc and the induced metric is fiat: 

jabdx°'dx'' — —Tcdr^ + dxidx^ 



(16) 



and hence satisfies the desired boundary condition. Here and henceforth x"" ^ {x^,t). The extrinsic 
curvature and unit normal on are 



1 



Tab labl'^'^Tcd 

P 



1 



dr 



1 - 



P 



dr + -=d, + 0{e^). (17) 



The Brown- York stress tensor is 
dxl 



Tnhdx^dx^ 



■ 2 o , ^I^^!l±IMdx-dx^ ^2^dx'dx^ +0{e^). 



dr^ -2—dx'dT 



(18) 



We wish to solve the Einstein equations as a power series in e. We first consider the necessary but 
not sufficient condition that the constraints be satisfied on Sc. At order €^ the metric is flat and Tab 
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is constant so they are trivially satisfied. The only way to get an order e term is with one power of 



and no derivatives. Such a linear term cannot appear because the constant terms in (14) can, through 



quadratic order, be obtained from a boost of flat space. The first nontrivial equation is encountered at 
order e^: 

This equation is satisfied if and only if is the velocity field of an incompressible fluid. Taking this to 
be the case, one finds at order e^: 

r^J^d'^Ta^ = drV, - 7jd^V, + d,P + d,v, = 0. (20) 

This is satisfied if and only if solves the Navier-Stokes equation with pressure P and viscosity 77 = r^- 
Once the constraints are satisfied it is ordinarily possible to evolve the solution off the hypersurface, 
in this case in the radial direction, at least for a finite distance. Here we have the danger of singularities 
at the horizon T-i^ near r = 0, or equivalently waves coming up from . We know from 28| that such 
singularities are absent in the linearized analysis provided the fluid viscosity takes the required value 



(12). We have checked by direct computation that this absence of singularities extends to the nonlinear 



case as well. That is all components obey 

Gra, Gab, Grr = C(e^) (21) 

and are nonsingular for finite values of r. Presumably the order and higher terms in the metric can 
be chosen so that the Einstein equations are solved exactly. 

It turns out that it is still possible to solve the Einstein equations analytically through order with 
the "wrong" value of the viscosity (i.e. rj ^ r^) even in the nonlinear case. As expected these solutions 
develop a singularity at r = near l-i^ , and are presented in Appendix [A] . 

5.2 Forcing the fluid 

Solutions of the linearized Navier-Stokes equation decay exponentially in the future. There is some 
expectation - although no proof - that nonlinear solutions eventually decay as well. Therefore the extrinsic 
curvature on in our examples is expected to become constant. 

On the other hand, already at the linear level, Navier-Stokes solutions grow exponentially in the far 
past and typically are singular at t = — oo . Therefore we expect that the dual geometry is also singular 
at r = —CXI, which is the past horizon "H^ of Sc. This singularity is not problematic for real fluids, as we 
are typically interested in cases where forcing terms correct the Navier-Stokes equation. For example we 
might consider a fluid which is initially at rest, stirred at time r = r* , and then left to evolve according 
to the unforced Navier-Stokes equation. 

In fact this kind of situation is also very natural to consider on the gravity side. Consider flat 
Minkowski spacetime with a flat metric and constant extrinsic curvature on the boundary Ec for r < . 
We then stir it at r = by momentarily perturbing the boundary condition that the induced metric 
on Ec be flat. This will send out a gravitational shock wave along t = and excite the geometry for 



T > . The result should be an appropriate gluing of ( 14 ) along a null hypersurface to flat Minkowski 
space. This is depicted in Figure [2j 

At the linear level, it is possible to explicitly construct the glued geometry describing this situation 
through order e^. The metric is 

ds^ = — rdr"^ + 2dTdr + dxidx^ 

- 2(1 - r/rc)vidx'-dT + (1 - r/rc) (5jW, + d^Vj) dx'^dx^ "'^{^ ^ d'^Vidx^dr 

-^(r-r,) (4:(l~r/rc)Fi + —ai]dx'dT~—l3,,dx'dx' +... (22) 
V rr. J rr. 
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Figure 2: On the Sc surface, prior to t = t*, all initial data is trivial. At r = t, , a gravitational shock v^ave arrives. The 
shock forces the fluid on Sc, and consequently the Vi is nontrivial on Sc after r*. 



where Fi is an arbitrary function of obeying diF^ 
functions of and related to Fi by 



0. Pij and ai (which is divergence free) are both 



djUi 



(23) 



Since the metric on Sc is no longer flat, the constraint equations become linearized Navier-Stokes with a 



forcing term similar to that described in 28 . For this configuration we have 



drv' - r^d^v' = F\x)5{t - n). 



(24) 



Clearly, since Wi(x,r) is taken to vanish for t < T^^,, the forcing term will cause it to jump to Fi{x^) 
at T = r*, after which it will evolve according to Navier-Stokes. Given (24 1 this geometry solves the 



linearized Einstein equations everywhere, and is characterized by an arbitrary divergence-free vector field 
Fi{x). Before r = r* it is flat, while afterward it is, up to a coordinate transformation, the linearization 



of (14) 



At the nonlinear level, the equations are cumbersome and we have been unable to explicitly construct 



the analog of (22 1 away from S^. However it seems plausible that qualitatively similar solutions persist 



at the nonlinear level. 



5.3 Singularities at r = ±00 



The square of the Riemann tensor for the solution ( 14 ) is given by 

7^2 = - A (^d^y^ _ d^y^f „ 2^ [d\,d^v' + 3dv {d.d^v^ - d,d\,)] + .... (25) 

This expression diverges at r = ±00. Of course perturbation theory cannot be trusted when \r\ is of order 
^, so the computation is unreliable in this regime. Whether or not there are actual divergences in these 
regions will depend on the details of the solution. In general, at r = — cx), black hole type singularities 
may plausibly arise. 

The divergence at r = -f-oo is outside the cutoff surface, so a priori need not concern us. Still we may 
ask what happens if we try to extend the solution to this region. In general relativity with no cosmological 
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constant it is hard to find solutions which are asymptotically flat in codimension one; i.e. there are no 
codimension one black holes. This suggests that many configurations will be singular if extended to 
r = +00. On the other hand, if we add a negative cosmological constant, there are codimension one 
asymptotically AdS black holes. At large r the cosmological term tends to dominate, and we expect in 
this case many solutions to have nonsingular extensions to this region. However, as we will see below, 
the hydrodynamic regime is small r so the large r behavior is of limited interest for the present purposes. 



5.4 Uniqueness 



Equation ( |14[ ) gives the first three orders in the e-expansion of metrics satisfying the Einstein equations 
with the prescribed boundary data. These solutions are constructed from nonlinear solutions of the 
incompressible Navier-Stokes equations. The latter are in turn, assuming existence and uniqueness for 
Navier-Stokes, specified by a divergence-free vector field ^'(a;, r*) at a moment of time r*. 

One may ask whether or not ( 14 1 is the unique solution with the prescribed boundary data (up to 
coordinate transformations and field redefinitions) associated to a given u*(x, r*). This can be addressed 
in the context of a combined weak-field expansion and e-expansion. The problem was solved to leading 
nontrivial order in the weak-field expansion in [28, . The unique solution is the first two lines of (14), but 
with a obeying the linearized Navier-Stokes equation. Generally one does not expect the dimension of 
the solution space in weak-field perturbation theory to change unless there is a linearization instability 
and associated obstruction. In the present case, the only potential obstruction is the Navier-Stokes 
equation which we are assuming can be solved. Hence one expects the solution ( 14 ) to be unique at each 
order in the e-expansion, up to the usual ambiguity of adding solutions of the leading order equations at 
subleading orders. 



6 Alternate presentation 



In this section we give an alternate presentation of the metric ( 14 ) in which all the factors of e appear 
explicitly, without being hidden in the functional dependence on the coordinates. This is accomplished 
by first transforming to hatted coordinates 



so that df = C(e") and we denote di 



TcX 

e 

a 



Hs'^ - 

"''p+2 — 



2Ti 



,1 - f 



~ di 

2rl 



T = — , r = TcT 

e^ 

©(e^). In the new coordinates 

^2 



(26) 



e^ e^ 



Vidx^dr — 2rrVidx^df 



+ {l-f) 



f 2P^.2 

dr + rcViVjdx dx-' 



+ rdv'^ + 2P)dTdf 



(f^ - lYcd'^Vidx^'dT + 



(27) 



(28) 



where P{x,t) = ■^P{x{x),T{f)), Vi{x,T) = ^Vi{x{x),T{T)), v' = iiiS'^^Vj and i,j indices are raised and 
lowered with Sij. The usual Navier-Stokes equation for u, P with t] — Vc implies 

dfVj - d'vj + v'^dkVj + djP = 0. 



(29) 



This is the Navier-Stokes equation with rj = 1 and no factors of e or 



Finally let us consider the rescaled metric ^5^+2 — ^d,^'p+2- The Einstein tensor is invariant under 



such constant metric rescalings. Rearranging terms and defining 



(30) 
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one finds 



+ [2dTdf + dx.dx' - 2(1 - f)vidx'dT + (1 - f){v'^ + 2P)df2] (31) 
+ A [(1 - r)viVjdx'dx^ - 2vidx'df + {v^ + 2P)dfdf + (1 - f^)d^v,dx'df] + .... 



The Brown- York stress tensor is 



1 



+ C'(A3/2). ^32) 

The important point here is that the geometry depends only on the ratio A = — and not e or Tc separately. 

Given that the rescaled geometry depends only on A and the e-dependence (21) of the unrescaled 
geometry ( 14 ) we conclude that in the hatted coordinates 

^2 „2 



C 

TT 7^T 



Gff 



2 ^rr 



0{X), Gff^rlGrr-0{X^) 



C = ^Gt 



(33) 



0(Ai/2), = ^G„ ^ 0(A3/2). 
go e 

Given the explicit factor of A~^ in gff, it is not immediately obvious in this presentation that in a direct 
computation the Einstein tensor will even have a good Taylor expansion in A. What happens is that, 
because gff = 0, there are only a limited number of powers of g-ff that can appear in the Einstein tensor, 
and one may thereby directly recover (33). In fact, direct computation reveals we do slightly better; the 
last line may be replaced by 



= -^Gt 



0(Ai), - ^G„ ^ 0(A2). 



(34) 



Notice that Gff in (33) is of order A'^ rather than A^. We can improve this by computing a few higher 
order pieces of the metric. Specifically, we add to (31) 



-2A(1 - f)q,dx'dT + 2X^gf^dfdx' + X^g^fdx'dx^ + ... 
Demanding that the r-independent pieces of Gff = be solved through order A*^ then fixes qi(T,x): 



(35) 



1 / - - \ " 



(36) 



Apparently qi is a kind of heat current. Demanding that the entire Gff — through order A" gives us a 
differential equation for the combination (5(f,f,x) = —2d^g^?^ + dfg^^'^^: 

Q+2fdfQ = 2d^q' -2v4 + 3fd^v^ - "^ (J),vj + djv^y + 2^^VJ^^v' -v^ djv^ + {^^^ (37) 

Choosing qi, Q accordingly, we find that as desired all components of the Einstein equations vanish for 
A -> 0: 

Gfa,G^-,,Gff^OiX). (38) 



7 Petrov type 

Interestingly, this geometry is of an algebraically special type. We consider the case of p = 2 to connect 
to the well-studied Petrov classification of spacetimes ||37j). A geometry is Petrov type II if there exists 
a real null vector fc'' such that the Weyl tensor satisfies 

W^,plAx]k''kP = 0. (39) 
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This happens if the invariant — 27 vanishes where I, J are both specific combinations of Weyl tensor 
components which can be found in ,37). For the metric (14), the lowest nonzero entries for /, J are at 
0{e'^) and O(e^) respectively. Hence, the first contribution to the invariant would be at 0{e^^); however 
the invariant vanishes through 0{e^^). At higher order in e, it gets modified by corrections to (14); we 
expect that including higher order terms in (14) enables (39) to be satisfied exactly. 



8 Nonlinear solution in the near-horizon expansion 

In section 4, the nonlinear Einstein equations with certain boundary conditions were solved in the non- 
relativistic, long-wavelength hydrodynamic e-expansion. This generalized the analysis given in 1281 of 



the e-expansion for linearized modes. 28 also considered, for linearized modes, a second, near-horizon 
expansion. Although physically inequivalent, the two expansions were found to be equivalent mathemat- 
ically and reduce to the linearized dynamics of an incompressible fluid. In this section, we consider the 
nonlinear version of the near-horizon expansion and find that it is again mathematically equivalent to 
the nonlinear e-expansion. 

In the e-expansion one solves the shooting problem for long-wavelength perturbations of Sc with a 
fixed leading-order extrinsic curvature. The proper acceleration of a worldline at fixed in Ec is to 
leading order just proportional to Ktt, so we may also view this as fixing the acceleration away from the 
origin. In the near-horizon expansion, instead of expanding in the wavelength one expands in the inverse 
acceleration. We begin with the flat metric on the Rindler wedge 

dsl_^^2 = -rdr'^ + idrdr + dxidx\ (40) 

To avoid confusion with the notation of the previous section we put the boundary on the accelerating 
surface 

r = fc, (41) 

so that r <fc. The near- horizon, large acceleration, limit is fc — 0. In order to exhibit the fc-dependence 
explicitly in the metric we transform to r = fcf, t = ^ so that r < 1 and 



ds1j^2 ^ -■^df'^ + 2dfdf + dxidx\ (42) 



In these coordinates the near- horizon limit rescales to infinity the coefficient of df^ at any finite f. 

We now wish to consider perturbations of this metric solving the Einstein equations order by order 
in the near- horizon expansion parameter that are consistent with a fiat induced metric at f = 1. At 
the level of linear perturbations, the most general solution was found in [28| (characterized in terms of 
the data at r = fc). This solution is (for all r) 



= - ^dr"^ + 2dfdf + dx^dx' - 2(1 - f)v^dx'df + [(1 - f^)d^Vidx'dT - 2vidx'df] + 0{fl) (43) 



where diV^ = and dfV^ — d'^v^ — 0. That is, is an incompressible fiuid flow obeying the linearized 
Navier-Stokes equation with unit kinematic viscosity. 



The nonlinear generalization of (43) which solves the nonlinear Einstein equations to 0{fc) is 



+ [2dfdr + dxdx' - 2(1 - f)v,dx'dT + (1 - f){v^ + 2P)df2] 

-I- fc[(l - f)viVjdx'dx' - 2vidx'-df + {v"^ + 2P)dfdf + (1 - f'^)d'^Vidx'dT - 2(1 - f)qi{f,f,x)dx'dT] 
+ rl (f, f, x)dx'dr + glf (f, f, x)dx'dx^] + Oirl) (44) 
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provided div'' — 0, dfVj 



(2) 
9fi ' 



glj are solutions of first order differential 



'vj + v'^dkVj + djP = 0. Qi 
equations of the type (36 1 and (37). Further 0{f'^) pieces do not affect the equations of motion to this 
order. 



We can now see explicitly that making the notation change w — >■ u, — )■ and f^, — )■ A in (44| 
gives us the rescaled solution (31) in section 6. Hence the near-horizon and hydrodynamic expansions 
are mathematically equivalent. 

Since we are identifying fc = A = Tc — > oo in the metric |l4[ ) is actually equivalent to the near- 
horizon limit f c — in ( 44 ) . This may at first seem odd, but the near-horizon- hydrodynamic equivalence 
involves a constant rescaling of ( 14 ) by a factor of -\ , the proper distance to the cutoff surface in the 



rescaled metric ( 44 ) indeed behaves as 
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A Appendix 

In the e-expansion. 



dSp^2 = ~ ''^dr'^ + Sdrdr -|- dxidx^ ~'^[^ ~ ) ^idx''dT — 2 — dx^dr 



+ 1 



{v^ + 2P)dT 



2 I "^dx'dx^ 



1 1? 


2P\ 







d^Vidx^dr — 2 ( 1 ) qidx^dr — 2ci (rlogr — TclogTc) d^Vidx^dr 



(45) 



2ci log ( ^ j {d,Vj + djv,) dx'dT + 2ci(l-^] v^djVidx'dr + ©(e"*) 



solves the Einstein equations through O(e^) if Vi obeys the incompressible Navier-Stokes equation with 
the "wrong" viscosity rj — Tc (1 + ci) where ci is a nonzero constant. For this geometry, the square of 
the Riemann tensor is 



ci (ci + 2) 



■ 1 2 

2diVjd^v'' + - {diVj - djVi) 



(46) 



which clearly diverges at r = unless ci vanishes or ci = —2. The last possibility is the time reverse of 
the first and exponentially growing in the future. 
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